The rate and luminosity function of Short GRBs 



o 

(N 



Tsvi Piran* and Dafne Guetta f 

*Racah Institute for Physics, The Hebrew University, Jerusalem 91904, Israe 
1 Osservatorio astronomico of Rome v. Frascati 33 00040 Monte Porzio Catone, Italy 

Abstract. We compare the luminosity function and rate inferred from the BATSE short hard bursts 
(SHBs) peak flux distribution with the redshift and luminosity distributions of SHBs observed by 
Swift/HETE II. The Swift/HETE II SHB sample is incompatible with SHB population that follows 
the star formation rate. However, it is compatible with a distribution of delay times after the SFR. 
D ' This would be the case if SHBs are associated with binary neutron star mergers. The implied SHB 

I^H [ rates that we find range from ~ 8 to ~ 30/iy Gpc _3 yr _1 . This rate is a much higher than what was 

ly-) ■ previously estimated and, when beaming is taken into account, it is comparable to the rate of neutron 

t— I ' star mergers estimated from statistics of binary pulsars. If GRBs are produced in mergers the implied 

rate practically guarantees detection by LIGO II and possibly even by LIGO I, if we are lucky. Our 
■ analysis, which is based on observed short hard burst is limited to bursts with luminosities above 

10 49 erg/sec. Weaker bursts may exist but if so they are hardly detected by BATSE or Swift and hence 
their rate is very weakly constrained by current observations. Thus the rate of mergers that lead to a 
detection of a gravitational radiation signal might be even higher. 
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^ INTRODUCTION 

The luminosity and rate of Gamma-Ray Bursts (GRBs) is one of the key issues in any 
astrophysical model. Shortly after the discovery that GRBs are cosmological Piran [1] 
and Mao and Paczynski [2] used the observed value of (V /V max ) = 0.32 and estimated, 
assuming that the GRB rate is independent of redshift and that they are standard candles, 
the rate of GRBs. Already at that time Piran [1] suggested that if GRBs arise due to 
neutron star mergers [3] their rate will depend on redshift and it will follow the neutron 
star binary formation rate with a logarithmic distribution of time delays, reflecting the 
distribution of merger times. 

In 1993 Koveliotou [4] noticed that the GRB distribution can be divided into two 
subsets of long and short bursts with a dividing duration of 2sec. As short bursts are 
also harder than long ones [5, 6], they are denoted as Short Hard Bursts (SHBs). Mao 
et al, [7] realized that BATSE is less sensitive to short GRBs than to long ones and 
pointed out that this should be considered when trying to fit the observed peak flux 
distribution to different models. Shortly afterwards Cohen and Piran carried out the first 
separate analysis of the luminosity and rates of long and short GRBs finding that the 
observed SHB population is located z < 0.5, namely they are much nearer than the 
observed distribution of long ones. Unfortunately a stubborn referee forced these authors 



1 talk given by Tsvi Piran 



to take out the discussion of short GRBs from [8] and these results were reported only 
in a conference proceedings [9] and in Cohen's Phd thesis [10]. These finding were 
corroborated shortly later by Katz and Canel [11] and later by Tavani [12] who found 
that (V /Vmax) of SHBs is significantly higher than the one measured for long bursts. 

With the discovery in 1997 of GRB afterglow and the subsequent identification of 
host galaxies and redshift measurements direct estimates of the luminosity and rates 
were obtained for long bursts. It was discovered that long GRBs follow the SFR and that 
typical (isotropic equivalent) peak luminosities are about 10 51 ergs/sec. 

However, until recently no afterglow was detected from any short burst who remained 
as mysterious as ever. Last spring Swift and then HETE II detected X-ray afterglow from 
several short bursts. In some cases optical and radio afterglows were detected as well. 
This has lead to identification of host galaxies and to redshift measurements. While the 
current sample is very small several features emerge. First, unlike long GRBs that take 
place only in star forming spiral galaxies SHBs take place also in elliptical galaxies 
in which the stellar population is older. In this SHBs behave like type la Supernovae. 
Second, the redshift and peak (isotropic equivalent) luminosity distributions of the five 
short bursts (see Table 1) confirm earlier expectations [9, 10, 11, 12, 13] that the observed 
SHB population is significantly nearer than the observed long burst population. 

We [13] (denoted hereafter GP05) have estimated, before the launch of Swift the 
luminosity function and formation rate of SHBs from the BATSE peak flux distribution. 
We have shown that the distribution is compatible with either a population of sources that 
follow the SFR (like long bursts) or with a population that lags after the SFR [1, 14]. 
If SHBs are linked to binary neutron star mergers [3] the SHB rate is given by the 
convolution of the star formation rate with the distribution P m {x) of the merging time 
delays x of the binary system. These delays reflect the time it takes to the system to 
merge due to emission of gravitational radiation. 
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Table 1: The Swift/HETE II current sample of SHBs with a known redshift. 

THE LUMINOSITY FUNCTION OF THE BATSE SHB SAMPLE 

As BATSE is less sensitive to short bursts than to long ones [7], even an intrinsic SHB 
distribution that follow the SFR gives rise to an observed distribution that is nearer to 
us. Still a delayed SFR distribution (that is intrinsically nearer) gives rise to even nearer 
observed distribution [13]. Therefore the recent observed redshift distribution of SHBs 
favors the delayed model and hence the merger scenario. Still the question was posed 
whether the predicted observed distribution is consistent with the current sample. Gal 
Yam et al., [15] suggested that the distributions are inconsistent and hence the delayed 
SFR model is ruled out. We re-examine the situation here and we show that while a 
delayed distribution with the best fit (maximal likelihood) parameters ia only marginally 
consistent (p K $ = 0.05 — 0.06) with the current sample, a delayed distribution with 



parameters within la from the best fit parameters is compatible (pks = 0.22 — 0.25). 
We discuss the implications of this result to the binary Neutron star merger model and 
to the detection of gravitational radiation from such mergers. We also discuss the recent 
suggestion of Nakar et al. [16] that there rate of mergers producing very weak bursts is 
very high. 

Our data set and methodology follow [17, 13]. We consider all the SHBs detected 
while the BATSE onboard trigger [18] was set for 5.5c over background in at least 
two detectors in the energy range 50-300keV. These constitute a group of 194 bursts. 
We assume the functional form of the rate of bursts (but not the amplitude). We then 
search for a best fit luminosity function. Using this luminosity function we calculate the 
expected distribution of observed redshifts and we compare it with the present data. We 
consider the following cosmological rates: 

• (i) A rate that follows the SFR (We do not expect that this reflects the rate of SHBs 
but we include this case for comparison.). 

• (ii) A rate that follows the NS-NS merger rate. This rate depends on the formation 
rate of NS binaries, that one can safely assume follows the SFR, and on the 
distribution of merging time delays, P m {x). This, in turn, depends on the distribution 
of initial orbital separation a between the two stars (t a 4 ) and on the distribution 
of initial eccentricities. Both are unknown. From the coalescence time distribution 
of six double neutron star binaries [19] it seems that i> m (log(T))<ilog(T) ~ const, 
implying P m (x) 1/t,[1]. Therefore, our best guess scenario is a SBH rate that 
follows the SFR with a logarithmic time delay distribution. 

• (iii) A rate that follows the SFR with a delay distribution P m (x)dx ~ const. 

• (iv) A constant rate (which is independent of redshift.). 

For the SFR we employ R S F2 of Porciani & Madau [20] : In models (ii) and (iii) the rate 
of SHBs is given by: 

ft(z) 
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Table 2: Best fit parameters Rate(z=0) , L*, a and /3 and their la confidence levels 
for models (i)-(iv). Also shown are the KS probability (p ks ) that the five bursts with a 
known redshift arise from this distribution. We show two results for KS tests one with 
GRB 0508132 at z = 0.7 and the other at z = 1.8. Case ii a corresponds to case ii with 
an L* value lower by 1 a than the best fit one. Other parameters have been best fitted for 
this fixed number. 



Following Schmidt [21] (see also [17, 13]) we consider a broken power law peak 
luminosity function with lower and upper limits, 1/Ai and A 2 , with power law indeces, 
a, /3 and luminosity break L*. 

We use Ai )2 = (30, 100) [13]. In [22] (denoted GP06 hereafter) we show that both 
limits are chosen in such a way that a very small fraction (less than 1%) of the observed 
bursts are outside these range. Hence one cannot infer anything from the observations on 
the luminosity function in this range. Comparing the predicted distribution with the one 
observed by BATSE we obtain, the best fit parameters of each model and their standard 
deviation. The results are shown in table 2 and in figs 1 and 2 of GP06. 

A COMPARISON WITH THE CURRENT SWIFT-HETE II SHB 

SAMPLE 

We can derive now the expected redshift distribution of the observed bursts' population 
in the different models. 

We assume that the minimal peak flux for detection for Swift is ~ 1 ph/cm 2 /sec 
like BATSE (note the different spectral windows of both detectors which makes Swift 
relatively less sensitive to short bursts). 

Fig. 1 depicts the expected observed integrated redshift distributions of SHBs in the 
different models. As expected, a distribution that follows the SFR, (i), is ruled out 
by a KS test with the current five bursts (pus < !%)■ This is not surprising as other 
indications, such as the association of some SHBs with elliptical galaxies suggest that 
SHBs are not associated with young stellar populations. 

A distribution that follows the SFR with a constant logarithmic delay distribution, (ii), 
is marginally consistent with the data (pks ~ 5 — 6%). The observed bursts are nearer 
(lower redshift) than expected from this distribution. If we use the Rowan-Robinson 
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FIGURE 1. A comparison between the expected integrated observed redshift distributions of SHBs for 
models (i)-(iv) and (ii a ) and the distribution of known redshifts of SHBs. 



SFR [13], rather than SF2 of Porciani-Madau, the situation is even more promising 
Pks ~ 10% (f° r either z=0.7 or z-1.8). When we move to a distribution that is la away 
from the best fit distribution we find p K $ ~ 22 — 25% and even higher for the RR SFR. 
Thus the suggestion [15] that the observed sample rules out the NS merger model (with 
a logarithmic merger time distribution) was somewhat premature. Note however, that the 
local rate with this model (ii a ) is sixteen times larger than the rate of the best fit model, 
(ii). This reflects the large flexibility in modeling the peak flux distribution. 

To demonstrate the flexibility of the data we have considered two other time delay 
distributions. A uniform time delay distribution (iii) and an overall constant SHB rate 
(iv). Both models are compatible with the BATSE SHB distribution and with the sample 
of SHBs with a known redshift {p K $ ~ 80% and 40% respectively.). This result is not 
surprising. The BATSE peak flux distribution depends on two unknown functions, the 
rate and the luminosity function. There is enough freedom to chose one function (the 
rate) and fit for the other. 

In all models compatible with the five bursts with a known redshift, the intrinsic 
SHB rate is pushed towards lower redshifts. The inferred present rates, ~ 30, ~ 8 and 
~ lO/z^Gpc^yr -1 for models (iii), (iv) and (ii<j) respectively, are larger by a factor 
ten to fifty than those estimates earlier (assuming that SBHs follow the SFR with 
a logarithmic delay with the best fit parameters GP05). The corresponding "typical" 
luminosities, L*, ranges from 0.1 to 0.7 x 10 51 erg/sec. 

CONCLUSIONS AND IMPLICATIONS 

We have repeated the analysis of fitting the BATSE SHB data to a model of the lumi- 
nosities and rates distributions. Our best fit logarithmic distribution model is similar to 
the best fit logarithmic model presented in GP05. A main new ingredient of this work 
is the fact that we consider several other models. We confirm our earlier finding that the 
BATSE data allows a lot of flexibility in the combination of the rates and luminosities. 

A second new ingredient of this work is the comparison of the best fit models to 
the small sample of five SwiftfHETE II SHBs. The SwiftfHETE II data gives a new 
constraint. This constraint favors a population of SHBs with a lower intrinsic luminosity 
and hence a nearer observed redshift distribution. It implies a significantly higher local 
SHB rate - a factor of ten to fifty higher than earlier estimates. The new observations 
of Swift show that the SHBs are nearer than what was expected before and therefore, 
their luminosity is lower and their local rate is higher. We stress that this new result was 
within the la error of the model presented in GP05, which had a very large range of 
allowed local rates and typical luminosities. 

Provided that the basic model is correct and we are not mislead by statistical (small 
numbers), observational (selection effects and threshold estimates) of intrinsic (two 
SHB population) factors we can proceed and compare the inferred SHB rate with the 
observationally inferred rate of NS-NS mergers in our galaxy [23] . This rate was recently 
reevaluated with the discovery of PSR J0373-3039 to be rather large as 80+^°/Myr. 
Although the estimate contains a fair amount of uncertainty [24]. If we assume that this 
rate is typical and that the number density of galaxies is ~ 10" 2 /Mpc 3 , we find a merger 
rate of 800+g9Q°/Gpc 3 /yr. Using a beaming factor of 30-50 for short bursts [25] this rate 



implies a total merger rate of ~ 240 — 1500/Gpc 3 /yr for the three cases (iii), (iv) and 
(iicr). The agreement between the completely different estimates is surprising and could 
be completely coincidental as both estimates are based on very few events. 

If correct these estimates are excellent news for gravitational radiation searches, for 
which neutron star mergers are prime targets. They imply that the recently updates high 
merger rate, that depends mostly on one object, PSR J0737+3039, is valid. These esti- 
mate implies one merger event within ~ 70Mpc per year and one merger accompanied 
with a SHB within ~ 230Mpc. These ranges are almost within the capability of LIGO I 
and certainly within the capability of LIGO EL If correct these estimates of the rate are 
excellent news for gravitational radiation searches, for which neutron star mergers are 
prime targets. 

To conclude we stress that we have assumed that the luminosity function has a lower 
limit of L*/30. This was just because even if such a limit does not exists weaker bursts 
would be barely detected. The current peak flux distribution of BATSE burst cannot 
confirm (or rule out) the existence of such population (note however, that Tanvir et al., 
[26] suggest that such a population exists on the basis of the angular distribution of 
BATSE SHBs). If such weak bursts exist then, of course, the overall merger rate will be 
much larger [16]. Such events will provide such a high rate that soon LIGO I will begin 
to constrain this possibility. 
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